Nonequilibrium Effects of Anisotropic Compression Applied to Vortex Lattices in 

Bose-Einstein Condensates 
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We have studied the dynamics of large vortex lattices in a dilute-gas Bose-Einstein condensate. 
While undisturbed lattices have a regular hexagonal structure, large-amplitude quadrupolar shape 
oscillations of the condensate are shown to induce a wealth of nonequilibrium lattice dynamics. 
When exciting an m = -2 mode, we observe shifting of lattice planes, changes of lattice structure, 
and sheet-like structures in which individual vortices appear to have merged. Excitation of an m = 
+2 mode dissolves the regular lattice, leading to randomly arranged but still strictly parallel vortex 
lines. 
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The experimental study of vorticity in a dilute-gas 
Bose-Einstein condensate (BEC) provides an useful per- 
spective on superfluidity. While the focus of early work 
was on the study of single or few vortices , recent ad- 
vances in technique have made it possible to create BECs 
containing large amounts of vorticity . In a rapidly 
rotating BEC, a lattice of vortices is the lowest energy 
state in the rotating frame 0. Experimental studies of 
lattices to date have concentrated on formation and decay 
processes [||-^,||J^. In this Letter we study the nonequi- 
librium behavior of vortex lattices under large-amplitude, 
anisotropic strain. 

Experimental technique For creating large amounts 
of vorticity in condensates, we use a modified version of 
the spin-up technique described in our previous paper 
l^. The starting point is a magnetically trapped cloud 
of 17 X 10'^ ^"^Rb atoms in the \F = l^mp = -1) state 
with a temperature approximately three times above the 
critical temperature for the BEC phase transition, Tc = 
67 nK. By manipulating the ellipticity of the confining 
potential, we can resonantly pump angular momentum 
into the as yet uncondensed cloud We then tune 

the trapping potential to near perfect axial symmetry 
with trap frequencies {ujp,uJz} = 27r{8.35, 5.45}IIz. In 
this round trap the cloud continues to rotate, and by 
performing a nearly one-dimensional evaporation along 
the axis of rotation (z-axis) as described in Ref. ||], we 
cool down and further spin up the cloud until a highly 
rotating condensate is born out of the rotating normal 
cloud. For the experiments described in this paper, the 
evaporation is continued to a point where little or no 
thermal fraction remains. As a result of this procedure, 
we obtain highly rotating condensates with typically 1 x 
10^ atoms, containing 130 vortices or more. 

To study these condensates we can image along the 
x, y, or z directions. When the condensate is held in the 
trap, the core size of the vortices is close to the resolution 
limit of our optical detection system so that individual 



vortices are not resolved by in-trap imaging. However, by 
nondestructively imaging a trapped condensate from the 
side (along the x-direction) and measuring its aspect ra- 
tio, the rotation rate of the condensate can be determined 
[^,||. For our highly rotating condensates we determine 
typical aspect ratios of 0.5, which is markedly different 
from the static aspect ratio of the magnetic trap, 1.53, 
because of centrifugal forces. Typical condensate rota- 
tion rates of 0.95wp are inferred. 

To resolve individual vortices, we can release a con- 
densate from the trap and image it after it has expanded 
by a factor of 5 [§,0. Such expansion pictures taken 
along the z-direction reveal a hexagonal lattice of vor- 
tex cores in highly rotating condensates [Fig. |l|(a)]. The 
lattice structure is remarkably regular even at the outer 
regions of the condensate, and is observed with very good 
contrast as shown by the cross section in Fig. 0(b). 

When looking at such vortex lattices in expansion from 
the side (i.e., along the x or y direction), good contrast 
is only obtained if the lattice is oriented in such a way 
that vortices line up behind each other along the direc- 
tion of view, or in other words when the direction of 
view is parallel to lattice planes. The vortex lattices are 
not stationary in the lab frame, but are rotating with 
the condensate, and furthermore the initial orientation 
of the lattice in the xy plane cannot be controlled in 
current experiments. Therefore the times when vortices 
line up along the direction of sight are unpredictable, 
and images revealing lattice planes are only obtained in 
a random subset of all trials. One example is shown in 
Fig. |l|(c), where the vortex planes are clearly visible as 
dark vertical lines. These images are direct verifications 
of the theoretical prediction that vortex lines are only 
marginally bent for our experimental parameters |]l3| , |l^ , 
as is also confirmed by the good contrast seen in the 
topview picture of Fig. 0(a). 

Vortex lattices such as the one shown in Fig. |l| are the 
starting point for our studies of the effects of dynami- 
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cally generated strain. We apply the distortion by reso- 
nantly driving quadrupolar shape modes in the vortex- 
filled condensate 16|. In the large-amplitude limit, 
the rriz = +2 and niz = —2 modes distort the circu- 
lar cross-section of the condensate into an ellipse and 
cause the major axis of the ellipse to rotate along with, 
or against, the sense of the condensate rotation, respec- 
tively. The generated strain is considerable - the ratio of 
major to minor axis of the condensate can be larger than 
three to one. The frequencies of the ruz = ±2 modes 
may be calculated using the sum rule argument given by 
Refs. Our condensates are rotating typically at 

O.OSwp, just 5% under the centrifugal limit. At the cen- 
trifugal limit, the major axis of the = +2 ellipse is 
fixed in the rotating condensate frame, and the — —2 
shape is fixed in the lab frame. At our rotation rate of 
0.95wp, the distortion, as seen in the frame of the lattice, 
rotates at only 0.4 Hz for the — +2 mode, compared 
with 8.34 Hz for the = —2 mode. For this reason 
we expect very different resulting lattice dynamics in the 
presence of either of the two modes. 

Lattice dynamics in the presence of an rriz — —2 
surface mode Because the = —2 mode is almost 
stationary in the lab frame for the parameters of our 
experiment, this mode can be excited conveniently and 
nearly resonantly by a static trap deformation. Indeed 
we observe that by jumping from a round trap to a trap 
that has an ellipticity of 3.6% ||2^, the ellipticity of the 
condensate increases over a timescale of 300 ms from ap- 
proximately to 40%, thereby exceeding the trap ellip- 
ticity by more than an order of magnitude. Since the 
vortex lattice is rotating quickly with respect to this al- 
most static deformation, it is an interesting and nontriv- 
ial question to ask whether or not a lattice structure is 
maintained, and, if so, how it rearranges in the presence 
of the excitation. 

Experimentally, vortex lattices can indeed be observed 
in expansion images up to 400 ms after the start of the 
continuously applied trap deformation (Fig. |^). Subse- 
quently, long range order is lost and vortex visibility be- 
comes low, presumably due to tilting []l6| or bending of 
the vortex lines jj]. A closer look at the first 400 ms of 
this evolution reveals a wealth of intriguing vortex dy- 
namics. As the vortex lattice rotates in the deformed 
condensate, the lattice planes must continuously shift rel- 
ative to each other to accommodate to the elliptical shape 
of the condensate. One consequence of this is transient 
changes of the lattice structure. In Fig. || the lattice has 
changed from the hexagonal structure of an undisturbed 
lattice to a near orthorhombic structure. 

The most striking observation, however, is of conden- 
sates containing sheet-like structures rather than individ- 
ual vortex cores, as shown in Fig. ^(a). We interpret these 
sheets as rows of vorticity along which individual vortex 
cores have essentially merged. Given the near perfect 
contrast of these sheets, as shown by the cross section in 



Fig. ||(b), we can exclude the possibility that the sheets 
are merely formed by a collective tilting of vortices along 
a lattice plane, which would result in an intermediate 
contrast. Similarly, we can rule out that this is an effect 
merely of imaging resolution - had the vortices retained 
their original structure but come so close together that 
we could not spatially resolve them, we again would see 
white stripes set off by grey troughs, rather than the near 
100% alternating stripes of white and black. If there con- 
tinue to exist wispy fingers of condensate that cross the 
stripes and differentiate individual vortices, they must be 
very tenuous indeed to be consistent with the observed 
contrast. The clouds have been continuously distorted 
from their original hexagonal symmetry, so the observa- 
tion of 16 stripes along the major axis of the BEG in 
Fig. ^a) tells us that there must also be the equivalent 
of 16 units of vorticity along the length of the stripes 
through the center of the cloud. By measuring the in- 
trap height of the condensate, the ellipticity in the xy- 
plane and the number of atoms, we can infer the ratio of 
unit vorticity spacing (along the minor axis) to healing 



length (given by 



where n is the mean density and 
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a the scattering length) to be on the order of 5. 

In fact the contrast of this sheet-like phase is so deep 
that the structure can be detected by looking at conden- 
sates in trap (along the x or z direction), even though 
the structure is close to our optical resolution hmit. By 
taking a time series of ten nondestructive images of a 
single condensate along the x direction, we see that the 
sheet structure appears and disappears periodically with 
a period of about 21.2 ms for a condensate rotating at a 
frequency of O.OSwp [Fig. 

This periodicity can be understood as follows: The 
sheet structure forms when a lattice vector lies along a 
minor axis of the cloud, presumably since this is when 
the vortices are closest together. As that vector rotates 
past the minor axis, structure along the vector is reestab- 
lished, the two-dimensional (2D) pattern of individual 
vortices reemerges, and the high-contrast stripes disap- 
pear from the in-trap images. A lattice vector aligns 
with the minor axis six times per lattice rotation period 
due to the sixfold symmetry of the unperturbed lattice, 
consistent with the observed frequency of periodic stripe 
formation. This interpretation is bolstered by the ex- 
pansion images, which show high contrast sheets when 
a lattice vector lies close to the minor axis, but 2D lat- 
tices (albeit distorted) for other orientations. We find it 
remarkable that this lattice-to-sheet-to-lattice sequence 
persists through multiple cycles before long range order 
and visibility are compromised. 

Using a trap with a stronger distortion of 34% as in 
Fig. ^c-e), we also observe sheets breaking apart in the 
center. In Fig. ^(c), the density distribution along the lat- 
tice vector that we see near-parallel with the minor axis 
is nearly featureless. Figure |^(d,e) hint at the succeeding 
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lattice evolution: Fig. ^(d) represents a small rotation 
from Fig. ||(c), and we see that the sheets have begun to 
reconnect across the cloud. Figure ||(e) appears to be a 
continuation of that process, with considerable evidence 
of uniformity in the repeating structure. 

Lattice dynamics in the presence of an = 
+2 surface mode We excite the — +2 mode by 
elliptically deforming the magnetic trap and rotating that 
ellipse around in the xy plane. When the uiz = +2 mode 
is driven approximately 1 Hz away from its resonance, 
the ellipticity of the condensate increases from to 50% 
over a timescale of 250 msec. Under the influence of an 
on-resonant drive, the ellipticity increases continuously, 
even passing beyond 90% after 700 msec and approaching 
100% As discussed above, the distortion rotates 

only very slowly with respect to the lattice, and thus we 
expect it to be correspondingly "gentle" with respect to 
its action on the lattice. 

Indeed, as we show in Fig. ||(a), we do not see the 
lattices squeezed into sheets for the niz — +2 mode. In- 
stead we observe a gradual increase of disorder in the 
lattice. 400 msec after we start the excitation, we see 
a complete loss of long-range order in the vortex struc- 
ture. The imaging contrast of individual vortices remains 
high, so we know they remain parallel to the line of sight 
(in contrast for instance to the case of Fig. ||), but oth- 
erwise they appear to have liquified (see the correlation 
histogram of Fig. This behavior confirms our expec- 
tation of the TTiz — +2 mode being a gentler disturbance 
for the lattice than the niz — —2 mode. 

Conclusions In conclusion we have succeeded in ob- 
serving a wealth of nonequilibrium vortex lattice dynam- 
ics in a BEC under anisotropic compression. Key ob- 
servations include the shifting of lattice planes relative 
to each other, leading to transient changes of the lat- 
tice structure, the formation of an unexpected and very 
pronounced sheet structure, and finally, a gradual tran- 
sition from an ordered lattice to a rather liquid-like ar- 
rangement in which vortex cores remain visible, but their 
relative positions are irregular. The observed sheet struc- 
ture is of particular interest: at present the behavior of 
densely packed vortex lattices is in the center of theo- 
retical discussions about vortices [^^-p^ . While most of 
these theoretical discussions propose to enter this regime 
by having high rotation rates almost inaccessible to ex- 
periment or inconveniently low atom numbers in the con- 
densate, this paper demonstrates an alternative and un- 
expected way of bringing vortices close together. We 
hope that the results presented in this Letter will stimu- 
late further theoretical discussions in the field of vortex 
lattices. 
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FIG. 3. Change of lattice structure, (a) Hexagonal struc- 
ture in an undisturbed lattice, (b) Near orthorhombic struc- 
ture seen transiently during lattice evolution in the presence 
of an mz = —2 quadrupolar surface mode. 
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FIG. 1. (a) Expansion picture of a vortex lattice seen along 
the rotation axis, (b) One pixel wide cross section along the 
white line in (a), (c) Expansion picture of a different conden- 
sate rotating more slowly than in (a), seen from the side. 
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FIG. 2. Lattice evolution after an rriz = —2 excitation. 
Pictures taken (a) 173 ms and (b) 873 ms after start of trap 
deformation 




FIG. 4. (a) Sheet-like structure as seen during lattice 
evolution in the presence of an = — 2 quadrupolar sur- 
face mode, (b) Cross section integrated over the white box 
in (a). Even though the box is wider than the calculated 
vortex core spacing, the observed contrast is nearly perfect, 
(c-e) Same as (a), but observed in a more deformed trap with 
= 27r{6.0, 8.6, 13.8}Hz 
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FIG. 5. (a-e) Nondestructive in-trap images of the 
sheet-like structures seen along the x-direction [conditions 
similar to Fig. ^(a)]. Spacing between images 10.6 ms. Note 
the very different spatial scale from expansion images (e.g.. 
Fig. |l]). (f) Cross section of (a), integrated over the conden- 
sate. 
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FIG. 6. (a) Irregular assembly of vortex cores as observed 
after 500 ms in the presence of an — +2 quadrupolar 
surface mode, (b) Histogram of distances between each pair 
of vortices in the regular hexagonal lattice shown in Fig. |l|a. 
The visibility of separate peaks reveals the high degree of 
long-range order, (c) Same as (b) but calculated for the BEG 
of Fig. Ma. Peaks are barely visible here, meaning that the 
lattice has almost reached a random disorder. 
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